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Abstract 


yThe  bounds  of  Boyd  (1959)  and  Soms  (1980a,  1980b)  for  the  tail  areas 
of  the  normal  and  t-distributions  are  used  to  obtain  a new  method  of 
evaluating  the  tail  areas.  The  absolute  and  relative  errors  and  numerical 
examples  are  given. 
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1 . Tht  Htthod 

Ht  btgin  by  Introducing  notation  and  stating  the  main  results  of  Boyd 
(1959)  and  Sons  (1980a,  1980b).  Boyd  (1959)  showed  that  If 

J 

♦(x)  ■ (2x)  ^exp(-x^/2),  T(x)  ■ ^(t)dt  , 


and  ■ T(x)/^(x),  x > 0,  then 


' "x  ‘ P<*-  l'--*  • 


max 


1 


where  p(x,v)  ■ (y  ♦ l)/[(x^  + (2/«)(y  + 1)^)^  * Yx],  y„.^  ■ 2/("-2), 

max 

Y„1n  ■ X - 1,  and  the  bounds  are  the  best  possible  In  the  class 
(p(x,y).Y  > - 1).  This  Is  also  discussed  In  Johnson  and  Kotz  (1970,  Ch.  33). 
Sons  (1980a,  1980b)  extended  the  above  results  and  showed  that  If  for  arbitrary 
real  k > 0 and  x > 0, 


Mt)  . c.mJsiUZil 


‘‘  r(k/2)(«k)’^^ 


Tk(x)  • >-F|,(x)  . f|((t)<lt, 

kk<*)  ■ 1^i(t«)/t(l+x*/k)f^(x)]  , 

'»'■  " » 2.  ■ <c^/(1-4c^)  .nd  Y„,„  ■ - 1 , and  for 


^ ■'min  '‘max  *'*®  Interchanged, 


and 


P(x,y) 


Ify 

(x^+4.c^(1+y)^)^^^+YX 


5C 


or  equivalently, 


(nf)fi,(x)p(x.Y^,„)  <7|^(x)  < • 

and  the  bounds  again  are  best  In  the  same  sense  as  for  the  normal. 

It  was  also  shown  there  that  If  k ■ 2,  ■ Y2  “ pCx.Yg). 

The  numerical  properties  of  these  bounds  are  discussed  In  the  above 
references.  The  Important  fact  to  be  noted  here  Is  that  the  bounds  control 
both  absolute  and  relative  error.  Using  the  bounds  as  a starting  point 
we  now  develop  a simple  method  of  evaluating  normal  ano  t>ta11  areas  that 
controls  both  absolute  and  relative  error,  as  opposed  to  the  usual  methods, 
which  generally  only  control  absolute  error. 

We  consider  estimates  of  the  tall  area  of  the  form 

for  the  tall  area  of  the  normal  and 

(?7K)P(»A„,„)f,(x)  . (1  - |i^)p(x,.„„)f,(x)  (1.2) 


for  the  tall  area  of  the  t.  We  want  the  estimates  to  He  between  the  upper 
and  lower  bounds  for  the  tail  area  and  be  strictly  decreasing  functions 
of  X and  therefore  Impose  the  added  restrictions  that 


be  > ad 
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and 


0 l|^37l  1 . all  X > 0. 


Since  f(0)  ■ we  may,  without  loss  of  generality,  assume  that  c » 1 and 
so  our  weight  functions  f are  of  the  type 


x/y\  . a»bx 
f(x)  - , 


(1.3) 


where0<a<_l,  d > 0,  be  > ad,  and  We  then  seek  that 

particular  choice  of  f which  minimizes  the  absolute  error.  A direct  computer 
search  led  to 


for  the  t,  where,  as  noted  before,  k Is  the  degrees  of  freedom.  (1.6) 

was  obtained  by  finding  the  optimal  constants  for  k = 25,  10,  5,  3,  1.5,  1,  .5 

and  fitting  a regression  line  to  them.  However,  In  the  interests  of  simplicity, 

for  k < 2,  we  did  not  Interchange  and  and  so  (1.5)  and  (1.6)  are 
— min  max 


understood  to  apply  for  all  k with  and  \ defined  as  for  k > 2. 
^ min  max 


Numerical  evidence  Indicates  that,  at  least  for  k = 1,  the  above  optimal 
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f(x) 

. .71 X 

"TT7T7 

(1.4) 

W" 

for  the  normal 

and 

• 

<W; 

f(x) 

b.x 

"W  • 

(1.5) 

'M 

^vv 

- .70  + 1.82/k  - .2/k^  , 

(1.6) 

b*  V • < 

•I*  ^ • 4 


estimate  Is  still  a decreasing  function  of  x. 


The  maximum  absolute  and  relative  errors  of  the  optimal  estimates  are 


remarkably  constant  over  the  range  1 1 k < ® and  hence  we  only  give  the 


normal  figures.  For  (1.4),  the  maximum  absolute  error  is  .66-10'^  and  the 


maximum  relative  error  is  .97-10"'^.  We  emphasize  once  more,  that,  unlike 


the  usual  methods,  which  generally  control  only  absolute  error,  the  above 


controls  both  absolute  and  relative  error  and  hence  can  be  used  to  calculate 


ordinary  and  Bonferroni  descriptive  levels  and  ordinary  and  Bonferroni 


percentiles. 


As  a check,  we  calculated  the  standard  textbook  table  of  the  normal. 


given,  e.g.,  in  Brown  and  Hollander  (1977)  and  found  at  most  a difference 


of  1 in  the  fourth  decimal  place.  We  also  compared  the  small  normal  percentiles 


given  in  Abramowitz  and  Stegun  (1965,  p.  977)  to  the  ones  obtained  from  (1.4) 


and  after  rounding  both  to  three  decimal  places  found  that  there  was  at  most 


a difference  of  1 in  the  third  decimal  place.  Similar  results  apply  to  the 


2.  Concluding  Remarks 


We  have  given  a method  of  calculating  normal  and  t-tail  areas  which  controls 


both  absolute  and  relative  errors.  The  listings  of  the  short  FORTRAN  programs 


are  available  on  request  from  the  author.  Preliminary  results  indicate  that 


it  is  possible  to  improve  on  the  accuracy  of  the  approximations  here  described 


at  a modest  increase  in  complexity  and  these  results  will  be  reported  shortly. 
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